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We present an analysis of the anomalous microlensing event, MOA-2010-BLG-073, announced 
by the Microlensing Observations in Astrophysics survey on 2010-03-18. This event was remark- 
able because the source was previously known to be photometrically variable. Analyzing the 
pre-event source lightcurve, we demonstrate that it is an irregular variable over time scales 
>200d. Its dereddened color, (V — I)s,o> is 1.221±0.051mag and from our lens model we derive 
a source radius of 14.7±1.3i?©, suggesting that it is a red giant star. We initially explored a 
number of purely microlensing models for the event but found a residual gradient in the data 
taken prior to and after the event. This is likely to be due to the variability of the source rather 
than part of the lensing event, so we incorporated a slope parameter in our model in order to 
derive the true parameters of the lensing system. We find that the lensing system has a mass 
ratio of g=0.0654±0.0006. The Einstein crossing time of the event, £e=44.3±0.1 d, was suffi- 
ciently long that the lightcurve exhibited parallax effects. In addition, the source trajectory 
relative to the large caustic structure allowed the orbital motion of the lens system to be de- 
tected. Combining the parallax with the Einstein radius, we were able to derive the distance to 
the lens, £>L=2.8±0.4kpc, and the masses of the lensing objects. The primary of the lens is an 
M-dwarf with M z , ] i=O.16±O.O3M while the companion has Mz, )2 =11.0±2.0Mj, putting it in 
the boundary zone between planets and brown dwarfs. 

Subject headings: Gravitational lensing: micro, microlensing, MOA-2010-BLG-073, MOA-2010-BLG- 
073L, brown dwarf, low mass binary 



1. Introduction 

The mass function of individual compact ob- 
jects (brown dwarfs and planets) in the Galaxy 
remains poorly understood, particularly at the 
low-mass end. Brown dwarfs (BD) are commonly 
defined as objects with masses between the deu- 
terium and hydrogen burning limits (DBL and 
HBL, respectively) but these can be hard to de- 
tect, being intrinsically faint and fading further 
as they cool over time. Below the DBL, the mass 
function for individual objects is even more poorly 
measured. Unbound, free-floating objects of plan- 
etary mass have been discovered via di rect imaging 



different period orbits. Brown dwarfs on the other 
hand, are considered to be the extreme low mass 
end of the star formation process by fragmenta- 
tion of locally o yer dense cores cause d by turbu- 
lence in a cloud ( Chabrier et al.ll201ll ) , which can 



them s elves form protopla netary disks ([Klein et al 
20031: IScholz et af1l2006l ). This mechanism may 



also produce object s of a few Jupit er masses. For 
a recent review, see iLuhmanl (|2012h . 



of clu sters (for example in a Orionis (jBejar et al 



20121 ) ) an d in the field (e.g. t h e 6-2 5 M, T object re- 
ported by iKirkpatrick et al ( 20061 )). Sumi et al 



( 201ll ) reported a population of planets which are 
either unbound or at very wide separations, dis- 
covered when their gravity caused short timescale 
microlensing events. 

At least partially as a result of these poor con- 
straints the origin of low-mass compact objects 
remain unclear. Although traditionally thought of 
as separate classes of objects, planets and brown 
dwarfs form a continuous scale of mass and are 
best disting uished by the circum s tances of their 



formation dBurrows et al 



20051 : ISahlmann et alJl2oTalchabrier et al.ll201lh 



2001: Chabrier et al 



Planets form in disks of material orbiting a pro- 
tostellar object and may subsequently migrate to 



The mass function of free-floating low-mass ob- 
jects i s likely to be different from those bound to 
stars. iMarcv fc But led (|2000l ) identified a paucity 
of BDs orbiting close (<3 AU) to their host stars, 
a region where Jovian-mass planets are commonly 
found. This "brown dwarf desert" may represent 
the gap between the largest objects that can form 
in protoplanetary disks and the smallest objects 
that can concurrently collapse / condense next to a 
star. 

Two different theories have been proposed to 
explain the forma tion of Jovi a n pla nets in disks 
(for a review, see IZhou et al.l (|2012ll V The core 
accretion model predicts planets form from pro- 
toplane tary cores, growing up to tens of Jupiter 
masses ( Mordasini et al. 20091 : Baraffe et al. .2010, 
e.g.) but predicts few giant planets around M- 
dwarfs. Higher mass stars are thought to have 
disks with enhanced surface densities which allow 



the cores to grow more rapidly ([Laughlin et al 
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2004), as do disks with a high fraction of dust, 
leading to enhanc ed planet format ion around high 
metallicity stars ( Ida fc Lin 2004 ). Alternatively, 
the model of planet for mation via g ravitational in- 
stabilities in the disk ([Boss 1120061 e.g.) tends to 



favor the formation of more massive planets, in 
generally wider orbits. 

A number of lines of evidence support the core 
accretion theory. There is a well-established cor- 
relation of incre asing planet frequency with stel- 



lar metallicity dSantos. Israelian fc Mavorl 2001 
iFischer k Valentil 120051 : iMaldonado et al.l \20m . 
The results of radial velocity surveys imply there 
is a dcrth of M-dwarf stars with massive, close- 
in planets. In part, this reflects an observational 
bias against these faint objects but the sample 
is sufficiently large that a real s t atistical trend 



is em erging (dimmin g et al.ll2008t Uohnson et al 



20101 : iBonfils et al.l 12011 ) for companions with 



P < 2000 d. Recent spectroscopic and Kepler 
results have confirmed the prediction of a rapid 
increase in frequency for planets with small radii 
(down to 2Rq) and P < 50 d for all spectral types, 
and found that these small planets are several 
times more common around stars of late sp ectral 
type (jBonfils et al.ll201ll iHoward et alll2012l ). 

However, the core accretion model has difficulty 
forming massive planets at large orbital radii and 
such syste ms have been disco vered, for example 
HR 8799 (jMarois et al.l 120081) . Furthermore, a 



number of planets have been found orbiting M- 
dwarf hosts at lar g er orb i tal separations, for ex - 
ample iDong et all (120061) : iForveille et all (|201lh : 



Batista et al.l (|2011l ). These systems may instead 



form through gravitational instability in the disk, 
which can account for companio ns up to sey eral 
Jupiter masses around M-dwarfs ( Boss 1 12006). 

To better understand the formation mecha- 
nisms of heavy substellar companions in bound 
systems we need to trace the distributions of the 
physical and orbital properties (such as mass ra- 
tio, orbital separation, occurrence frequency) of 
a significant number of systems. Yet relatively 
few bound brown dwarf companions have been re- 
ported, despite their being easy to detect at close 
orbital separations (the "brown dwarf desert"). 

Microlcnsing offers a complementary window 
onto BD and planet formation by probing for 
cooler companions of all masses in orbital radii be- 
tween ^0.2 - 10 AU, separations which are difficult 



or time consuming to explore by other methods 
( Shin et al.ll2012al ). It can probe the companion 
mass function down to M- and brown dwarf hosts, 
and is sensitive to companions from nearly equal 
mass down to terrestrial masses. 

Sixteen systems have been published to datc0, 
and thanks to large-scale galactic lensing sur- 
veys and efficient follow-up, each season's Bulge 
observing campaign is now pro ducing a regu- 



lar yi el d of new d i scover i es (e.g. iBachelet et al.l 



(l2012h : lYee et al.l (l2012h : iMivake et all (<2012h ) 



Of these 16 companions, 3 are giant planets or- 
biting M-dwarf s tars: OGLE-200 5-BLG-071Lb, 
a 3.8 Mj planet (|Dong et al.l 12009^ . MOA-2009- 
BLG -0387Lb, with M P =2.5 Mj planet (jBatista et al 
201ll ) and MO A-20 1 l-BLG-293Lb which hosts a 
2.4 Mj companion (|Yee et al.ll2012f ). 

Here we present the newly discovered system, 
MOA-2010-BLG-073L, an M-dwarf star with a 
companion whose mass is close to the deuterium 
burning limit of ~42.6 Mj. The discovery and 
follow-up observations are described in §[2] and we 
discuss the variability of the source star in § [3] We 
describe our analysis in § 21 from which we derive 
the physical properties of the lens in § [SJ Finally, 
we discuss our findings in § [6] 

2. Observations 



The microlensing event, MOA-2010-BLG-073, 
was first announced by the Microl ensing Ob- 
servat i ons in Astr o physi cjl (MOA, iBond et al 



(|200ll ): ISumi et al.l (|2003l ) on the 1.8 m telescope 
at Mt. John Observatory, New Zealand) survey 
on 2010-03-18. A background source star in the 
Galactic Bulge, a=18:10:11.342, £=-26:31:22.544 
(J2000.0), previously having a mean baseline mag- 
nitude of / ^16.5 mag, was discovered to be rising 
smoothly in brightness consistent with a point- 
source, point-lens (PSPL) microlensing event. 
However, on 2010-05-03 the event was found to 
show an anomalous brightening of ~ 0.5 mag and 
an alert was issued (K. Furusawa, private comm.). 

Microlensing follow-up teams worldwide 



listed on exoplanet.eu 

2 www. phy s . canterbury, ac . nz/ moa 
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RoboNet-lH dTsapras et al.ll2009h. uFUN 1 ! d Gould et al 



2006), PLA NET] 5 ! dBeaulieu et al 



MiNDSTErll (jDominik et all l2010( ) 



1200611 and 
- responded 

to provide intensive coverage of the event for the 
duration of the anomaly (^2 days), and monitored 
the event as it returned to baseline, over the course 
of the next ^2 months. 

In addition to the MOA data, taken with the 
wide band "MOA-Red" filter (corresponding to 
R+I bandpasses), the event was observed from 
several other sites in New Zealand. The 0.41 m 
telescope at Auckland Observatory, the 0.35 m at 
Kumeu Observatory and the 0.41m Possum Ob- 
servatory all used i?-band filters while the 0.304 m 
Molehill Astronomical Observatory (MAO), the 
0.35 m telescope at Farm Cove Observatory (FCO) 
and the 0.4 m telescope at Vintage Lane Observa- 
tory (VLO) all observed it unfiltered. The event 
was then picked up from three sites in Australia, 
firstly in /-band by the 1 m Canopus telescope in 
Tasmania followed by the 2 m Faulkes Telescope 
South (FTS), where an SDSS-i filter was used. 
The 0.6 m telescope in Perth also observed in / 
band. Of the observing sites around longitude 
zero, the event was imaged from the 1 m telescope 
at the South African Astronomical Observatory 
(SAAO) using an I band filter and in J, H and 
K s by the 1.4 m Infra Red Survey Facility (IRSF), 
also at SAAO. The 2 m Liverpool Telescope (LT) 
observed the event in SDSS-? from the Canary Is- 
lands. 

As darkness fell in the Americas, a number of 
Chilean telescopes picked up the observing baton: 
the SMARTS 1.3 m at the Cerro Tololo Interamer- 
ican Observatory (CTIO) obtained data in V, I 
and H bands with the ANDICAM camera, and the 
Danish 1.54 m used an / filter. Though in the 
midst of commissioning the new OGLE-IV cam- 
era at the time, the 1.3 m Warsaw telescope also 
covered the event in J-barufl. The 1 m Mt. Lem- 
mon Telescope in Arizona imaged the event in the 
/-band and in the extreme west, the 2 m Faulkes 
Telescope North (FTN), Hawai'i used an SDSS-i 



3 robonct. lcogt.net 

4 www. astronomy.ohio-state.edu/~microfun 
5 planet. iap.fr 

6 www. mindstep-science.org 

7 ogle . astrouw .edu.pl 



filter to complete the 24-hour coverage of the event 
from the Pacific. Table [T] summarizes the data ob- 
tained, which are plotted in Figure [T] 

The high density of Galactic Bulge star fields 
and the consequent degree of overlap (or blending) 
in stellar point-spread functions (PSF) has long 
since made difference image analysis (DIA) the 
photometry method of choice among microlens- 
ing teams. Both MOA and OGLE make their 
photometry available to the community, auto- 
matically reducing their data with t heir custom 



pipelin es d escribed resp e ctively in iBond et al 



(|200ll ) and lUdalski et all (|2003l ). The RoboNet 
data (from FTN, FTS and the LT) were re- 
duced with the project's automated data reduc- 
tion pipe line, which is b ased around the DanDIA 
package (jBramich I l2008h . This software was also 
later used to reduce data from Canopus, the Perth 
0.6 m, the SAAO 1 m and the Zf-band data from 
CTIO, while the DIAPL package was used to pro- 
cess the images from the Danish telescope. The 
PLANET team released their photometry (pro- 
duced by the WISIS pipeline) in real time via their 



webs ite, and the Pysis DIA pipeline (jAlbrow et al 



20091) was used for later re-reduction of these data 
sets. 

3. Variability of the Source Star 

MOA-2010-BLG-073 was present in the fields 
of the OGLE-II and OGLE-III surveys so the 
source star's I-band photometric record extends 
from 1998 to 2006 (Fig. [2J. OGLE-IV was in 
the commissioning phase when this event took 
place. From this excellent baseline it was immedi- 
ately clear that the source is variable over many- 
month time scales. This raised the possibility that 
shorter-term variability might obfuscate the mi- 
crolensing signal, making it difficult to determine 
its properties. 

To investigate this possibility, we performed a 
search for periodicities in the baseline OGLE-II 
and OGLE-III data, excludi ng the lensing event, 
using the ANOVA algorithm (jSchwarzenberg-Czernv 
1996(). Due to the seasonal gaps in the baseline, 



we analyzed the OGLE-II data in yearly subsets 
as it is the best sampled, searching for periods 
between P=0. 5-200 d. As Figure demonstrates, 
there are no significant or persistent periodicities, 
other than the expected integer multiples of the 
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Table 1 
Summary of observations. 



Telescope & 


Filter 




ry 


N frames 


N frames 


(<0 

a 


(c) 


j- r i 
aperture [mj 








total 


used 




r i 
[mag] 


MOA 1.8 


R/I 


0.7027 


0.6118 


1747 


1726 


1.305 


0.005 


OGLE 1.3° 


I 


0.6098 


0.5103 


47 


42 


2.188 


0.005 


Auckland 0.41 


R 


0.7027 


0.6118 


136 


136 


0.910 


0.000 


Canopus 1.0 


I 


0.6098 


0.5103 


162 


159 


1.310 


0.005 


CTIO 1.3 


V 


0.7817 


0.7048 


19 


18 


0.603 


0.000 


CTIO 1.3 


I 


0.6098 


0.5103 


162 


162 


1.010 


0.003 


CTIO 1.3 


H 


0.4145 


0.3206 


586 


575 


1.340 


0.016 


Danish 1.54 


I 


0.6098 


0.5103 


498 


491 


1.130 


0.014 


Farm Cove 0A b 


Unfiltered 




0.5611 


225 


225 


0.975 


0.000 


FTN 2.0 


SDSS-i 


0.6098 


0.5103 


159 


158 


1.055 


0.011 


FTS 2.0 


SDSS-i 


0.6098 


0.5103 


129 


129 


1.125 


0.006 


IRSF 1.4 


J 


0.4836 


0.3844 


4 


4 


1.000 


0.000 


IRSF 1.4 


H 


0.4145 


0.3206 


4 


4 


1.000 


0.000 


IRSF 1.4 


K S 


0.3550 


0.2684 


4 


4 


1.000 


0.000 


Kumeu 0.35 


R 


0.7027 


0.6118 


272 


272 


0.772 


0.000 


Lemmon 1.0 


I 


0.6098 


0.5103 


116 


105 


1.290 


0.020 


LT 2.0 


SDSS-i 


0.6098 


0.5103 


167 


167 


1.155 


0.006 


MAO 0.304 6 


Unfiltered 




0.5611 


238 


238 


1.025 


0.018 


Perth 0.6 


I 


0.6098 


0.5103 


66 


66 


1.095 


0.008 


Possum 


R 


0.7027 


0.6118 


15 


15 


1.030 


0.009 


SAAO 1.0 


I 


0.6098 


0.5103 


30 


30 


1.050 


0.011 


Vintage Lane 0.4 6 


Unfiltered 




0.5611 


124 


124 


0.995 


0.000 



includes only OGLE-IV data taken during event MOA-2010-BLG-073. 

b For unfiltered or very broadband observations we adopted a limb-darkening parameter 
value which was the average of that for R and I bands: (r_R + r_f)/2. 

c u\, T\ are defined by Eqns.[3l[4][5]and ao, ai by Eqn.(6] 
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Fig. 1. — Top panel: The complete light curve of MOA-2010-BLG-073 with the data sets from the various 
observatories overlaid with our best fitting model. The inset shows the anomaly in greater detail. Lower 
panels: The residuals of the fit for each of the best-fitting models of each class. 
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15.24 




500 1000 1500 2000 2500 3000 3500 4000 
HJD-2450000.0 

Fig. 2.— All available OGLE-II (1998-2001) and OGLE-III (June 2001-2006) lightcurve data for the source 
star of MOA-2010-BLG-073, taken prior to the event. The instrument upgrade to OGLE-III occurred around 
HJD=2452000.0. The apparent offset in magnitude at this time may be due in part to the difficulties of 
accurately calibrating photometry between the two wide field surveys. 




Fig. 3. — Periodograms of the OGLE-II (1998-2001) seasons of data, analyzed separately with the ANOVA 
algorithm. No significant periodicities were found in the range P=0. 5-200 d, aside from the expected aliases. 
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1-day sampling alias. 

We then combined the OGLE-II and III data 
sets in order to search for periods up to P=4000d. 
Figure [5] indicates a slight (~0.04mag) magnitude 
offset between the OGLE-II and III data. This 
can occur as a residual of OGLE's photometric 
calibration between the two surveys but it might 
also be the result of the intrinsic stellar variation. 
Therefore, we performed a search for periods be- 
tween 0.5 and 4000 d based on the combined data 
both with and without this offset (estimated vi- 
sually) . In both cases the periodogram was domi- 
nated by the window function; the only significant 
power was found in the peaks marking multiples 
of the 1-day alias, plus one peak at extreme low 
frequency corresponding to the finite length of the 
data set. We conclude that this star is an irregu- 
lar, long-term variable, most likely as a result of 
pulsations. 

However, there remained the possibility that 
the star could be irregularly variable on time scales 
comparable to that of the lensing event. To test 
this possibility, we binned the OGLE-II light curve 
on a range of time scales between 2-1200 d. To 
each bin we fitted two functions: one of constant 
brightness, and one with a linear slope and calcu- 
lated the weighted root-mean-square (RMS) pho- 
tometric scatter around each function per bin. We 
plot the average and maximum RMS (calculated 
over all bins in the lightcurve) against the width of 
the bins in time in Figure|3J On time scales shorter 
than 200 d the RMS scatter in the binned light 
curve is reasonably constant implying no signifi- 
cant short term variability. The longer term trend 
becomes clearly evident in the constant bright- 
ness curves for time scales longer than 400 d. We 
note that for bin widths between ~800-1200d, 
this curve has an RMS actually exceeding that of 
the whole light curve; this is because these bins 
were sufficiently wide that the first bin included 
the majority of the data, and the most variable 
sections of the light curve. As the bin widths be- 
came longer, they included more data points from 
the relatively stable section towards the end of the 
OGLE-II data set, and the RMS drops. The de- 
viation of the OGLE lightcurve from a constant 
brightness exceeds 3a for timescales longer than 
~750d. However, it deviates from a linear slope 
by < 1.6cr for timescales less than 1000 d, so we rep- 
resent this variation as a gradient in the lightcurve 



over the duration of the event. 

The most intuitive way to account for the vari- 
ation of the source was to measure the gradient of 
the lightcurve taken at baseline immediately be- 
fore and after the event. Unfortuantely only one of 
the available datasets covered these periods. Fit- 
ting a straight-line model (via a non-linear least- 
squares Marquardt-Levenberg algorithm) to the 
MOA 2009 and 2011 season data, we measured a 
slope of 0.018 mag/yr. However, the RMS scatter 
in the residuals of this fit were 0.016 mag, making 
it difficult to properly determine the slope. Ad- 
ditionally, in order to remove this slope from the 
other datasets, it would be necessary to transform 
the fluxes measured by each telescope on to the 
same scale as the MOA data. We attempted this 
via a linear regression approach but found that sig- 
nificant residuals remained. These contributed to 
overall higher y 2 values when the corrected data 
were fit with binary lensing models. We therefore 
adopted the alternative method of incorporating 
the slope as an additional parameter in our lens- 
ing model which we fit to the original, uncorrected 
data and we describe this approach in the follow- 
ing sections. 

4. Analysis 

In the analysis of this event, we used the es- 
tablished mo deling software developed by S. Don g 
and C. Han (|Dong et al.ll2006l : IShin et al.ll2012bh . 



4.1. Initial Parameters 

As a starting point for our analysis, we needed 
approximate values for the three parameters of 
the standard model for a PSPL event (not yet in- 
cluding the slope; this is discussed in § 14. 8j) : to, 
the time of peak magnification occurring at the 
closest projected separation between the lens and 
source, uq and the Einstein radius crossing time, 
tE- Following standard convention, all distances 
are quoted in units of the angular Einstein radius, 
0b, of the lens. 

To estimate these parameters, we combined all 
available data sets into a single light curve using 
the following scaling to take account of the vary- 
ing degrees of PSF blending from different instru- 
ments: 



f(t,k) = A(t)f s (k) + f b (k), 



(1) 
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where f(t, k) is the measured flux of the target at 
time t from data set k, A(t) is the lensing mag- 
nification at that time, f a (k) is the flux of the 
source star and /&(&) represents the flux of all stars 
blended with the source in the data set. A regres- 
sion fit was used to measure f s and ff, for each data 
set, producing an aligned light curve. Although 
the resulting parameter estimates are somewhat 
different from their 'true' values due to the ex- 
istence of the anomalous deviation on the light 
curve, they provided a starting point in parame- 
ter space. 

We note that two additional parameters can 
contribute to a PSPL model. They are the lens 
parallax parameters 7Te,jv and tte,e that account 
for the light curve deviation caused by the mo- 
tion of the Earth in its orbit over the course of 
the event. The vector microlens parallax, 7Te = 
AU/rE, where rg is the Einstein radius projected 
onto the observer plane and: 



(te.jVjTTe.b) = (cos0 w ,sin^)7r E , (2) 



where <f>^ represents the direction of lens motion 
relative to the source as a counter-clockwise angle, 
north through east. However, this is generally sig- 
nificant only for long (^months) time scale events, 
and was included at a later stage (see Section |4~6"|) . 

4.2. Finite Source 

The sharp spike feature in the lightcurve is in- 
dicative of the source closely approaching or cross- 
ing a caustic. In these circumstances, it cannot be 
approximated by a point light source and must be 
treated as a disk of finite angular radius, p, with 
wavelength-dependent limb-darkening. This is ad- 
dressed wi thin our software us ing the ray-shooting 
approach ( Kavser et al. 19861) : the path of light 



rays is traced from the image plane back to the 
source, taking into account the bending of the tra- 
jectory according to the lens equation. If a ray is 
found to "land" within the radius of the source, 
its intensity is computed taking limb-darkening 
into account. We derived this from the linear limb 
darkening law: 

7 A (cos0) = 7 X (1) [1 - u A (l - cos0)] , (3) 
where I\ is the intensity of the source at radius 



<fi from the center, relative to the central intensity 
in the same wavelength, A, scaled by the co- 
efficient u\. While more accurate limb darkening 
models are available, they are not commonly used 
in microlensing analyses due to the complexity in- 
troduced by combinin g data from many sources 
(|Bachelet et al.l f|2012h discussed this in more de- 
tail). The values of u\ for each passband were 
calculated from the Kurucz ATLAS9 stellar at- 

(Il979h us- 



mosphere models p resented by 



ing the method of iHevrovskv 



Kurucz 



(2007). However, 



within the microlensing community and software, 
Equation [3] m ore commonly fol lows the formalism 
derived bv TAlbrow et al.l (| 19991 ) : 



h = 



Fx_ 

7T6I2 



1 



1 



COS( 



(4) 



where F\ is the total flux from the source in 
a given passband and <j> is the angle between the 
line of sight to the observer and the normal to the 
stellar surface. The limb darkening coefficient, T\ 
is related to u\ by: 



2u> 



(5) 



The values of u\ and T\ applied for each 
dataset are presented in Table Q] The lensing 
magnification is then computed as the ratio of the 
number of rays reaching the source plane relative 
to the number in the image plane. This approach 
is only required while the source is close to the 
caustic. At larger separations, the software em- 
ploys a semi-analytic hexadecapole approximation 
to the finite so urce calculation to improve compu 
tation speeds ( Peicha fc Hevrovskv 20091 : Gould 
20081) . 



4.3. Standard Binary Model Grid Search 

To model the light curve of a binary lens event, 
we introduced three additional parameters: q = 
Ml,2/Ml,i, the ratio of the masses of the two 
bodies composing the lens where Mli is the more 
massive component, Sq, the projected separation 
of those masses and a , the angle of the trajectory 
of the lensed source star, relative to the lens' bi- 
nary axis. The frame of reference was defined to 
be at rest with respect to the Earth at time io,par, 
which we took to be the time of caustic crossing at 
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HJD=2455321.0, estimated from the easily iden- 
tifiable fea ture in the light cur ve (following the 



notation of S kowron et al.l l2011) 



With seven variables in the model (to, Uo, £e, 
ccqj s 0j Q, P)i a number of different lens/source con- 
figurations may produce similar light curves, so it 
was necessary to thoroughly explore a large area 
of parameter space in order to ensure all possible 
solutions are identified. We therefore constructed 
a grid of models, spanning set ranges in the val- 
ues of the three variables upon which the overall 
X 2 of the fit depended most sensitively, sq, q and 
ao- Each node in this grid took fixed values of 
(so,q,ao) and used a Markov Chain M onte Carlo 
(MCMC) approach (|Dong et al.ll2006l) to find the 
best fitting model by optimizing the other param- 
eters. To improve efficiency, a magnification map 
is generated by rayshooting for each point in the 
grid from which the model light curves used to 
compute the x 2 are drawn. The grid covered the 
following range: log(so) = —0.6 : 0.6 in steps of 
0.012, log(q) = -4.0 : 1.0 in steps of 0.05 and 
a = 0.0 : 6.3 in steps of 0.6. 

Mapping out the x 2 f° r each node in this grid, 
we found a number of local minima. Visual inspec- 
tion of these models overlaid on the light curve 
demonstrated that some more closely followed the 
data than others. Our first pass analysis included 
substantial baseline photometry before and after 
the event. This was not well fit by the models due 
to the variability of the source and hence the \ 2 
map gave a distorted view of regions in param- 
eter space that best match the event. For this 
reason, we proceeded by repeating the grid search 
using just data taken during 2010. This produced 
two clear minima in \ 2 of which one model stood 
out as by far the best match to the data. We 
then conducted a refined grid search over this re- 
stricted region of parameter space, taking smaller 
incremental steps. 

4.4. Optimized Standard Binary Model 

The refined grid search produced a reasonable 
model, fitting the majority of the data from all 
telescopes. This was used as a guide to identify- 
ing likely outlying data points, for which the qual- 
ity of the reduction was then checked. A handful 
of data points were removed at this stage. How- 
ever, this model included only fixed values for So, 
q and ao. To properly determine the standard bi- 



nary model for this event, our next step was to 
allow the seven parameters (io, uoi ^e, «o, s o, 
q, p) to be optimized during the MCMC fitting 
process, which used the grid search results as its 
starting point. At this point we included the ex- 
tended baseline data from MOA for seasons 2009 
and 2011, as these fall within the period for which 
the source's variability can be approximated with 
a straight line; we address this in Section 14.81 

4.5. Normalization of Photometric Errors 

When fitting microlensing events, the reduced 
X 2 of the fit on a per data set basis, X 2 e d> typ- 
ically produces a range of values both less than 
and exceeding the expected unity value. This can 
occur as different groups have slightly different 
ways of estimating photometric errors, but can 
lead to over- or under-emphasis being placed on 
particular data sets during the modeling process. 

A common technique to address this issue is to 
arrive at a complete model for the event and then 
use this model to renormalize the original photo- 
metric errors of each data set, e or i g , according to 
the expression: 



ao 



(6) 



We first conducted the sequence of models de- 
scribed in the following sections in order to find 
the best model for the event. We then set the 
coefficients ao, ai such that the Xrcd relative to 
that model equaled unity; the adopted values are 
given in Table [T] We then repeated our MCMC 
fitting process, starting with the Standard binary 
model and systematically adding parameters in to 
determine the extent of improvement in the model 
X 2 in each case. We compare models in Tables [5] 
& 131 and we plot the residuals (data — model) in 
Figure [TJ In the following sections, the x 2 values 
given are those post-renormalization. 



4.6. Parallax 

Given that the event's ^E~44d~0.12yr, it was 
necessary to include parallax in our model. Us- 
ing the parameters of the Standard Binary Model 
as a starting point, we allowed our fitting process 
to optimize for 7Te,_e, tte,n also. We found that 
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this significantly reduced the \ 2 °f the overall fit 
to 6064.773. By default, this procedure explored 
models with positive projected separations at clos- 
est approach of the lens and source, that is uq > 0, 
which we define as the source's trajectory passing 
the caustic at positive values of 9 y in the lens plane 
(see Figure [6]). For the standard model case the 
symmetry with respect to the binary axis of the 
caustic means that the uq < solutions are iden- 
tical. However once parallax was included, this 
was no longer the case, so we also explored uq < 
solutions (this de generacy is further discussed in 
Park et al The parameters of the uo > 

model were taken as a starting point for the fit, 
except that the sign of uq was reversed and the 
ao value became 2ir — a n . We found the best fit- 
ting u < model to be slightly less favored, with 
X 2 =6099.487. 

4.7. Lens Orbital Motion 

The mass ratio and projected separation de- 
termined from this model put this event close 
to the boundary between close and intermedi- 
ate/resonant caustic structure. In this regime, 
small changes in the projected separation of the 
lensing bodies due to their orbital motion can 
effectively change the shape of the caustic (see 
Fig. [6]) while the event is underway, sometimes 
causing detectable deviations in the lightcurve. To 
explore this possibility, we included additional pa- 
rameters in our model to describe the change in 
projected binary separation, ds / dt and the rate of 
change of the angle of the projected binary axis, 
da/dt. Again we found that the uq > model 
gave the best fit, with x 2 =5544.071, compared 
with 5690.003 when m < 0. This type of orbital 
motion is classi fied as "separational" in the schem a 
put forward bv lGaudi I ([20091 ) ; IPennv et al.l (|2010l ) , 
and is detected in this event as the source happens 
to cross the cusp of the caustic in the position 
where the caustic changes most rapidly. 

4.8. Sloping Baseline 

While taking these second-order effects into ac- 
count significantly improved the fit to the data, 
the overall x 2 remained rather high. Visual in- 
spection of the light curve still showed a gradi- 
ent, especially in the 2010 baseline before, rela- 
tive to after, the event. Based on our analysis in 



§ G3 this trend is likely to be part of the longer- 
term variability of the source and not associated 
with the microlens. In order to determine the true 
lens/source characteristics though, this trend had 
to be taken into account. 

The OGLE-II and III data demonstrate that 
the source variation over the *~150 d time scale of 
the event can be approximated by a straight-line, 
rather than a higher-order function. We there- 
fore introduced a 'slope' parameter to our model, 
representing the linear rate of change in magni- 
tude during the event. This further improved the 
X 2 , and the best-fitting model was once again the 
u > solution with x 2 =4782.367. 

We note that there exist degeneracies between 
the slope parameter and those for parallax and 
orbital motion as they can be used to fit simi- 
lar residuals in the lightcurve. To test for pos- 
sible degeneracies, we also fit a standard model 
plus the slope parameter alone and found that 
X 2 =5334.378. The value for the slope from this 
model, -0.0160±0.0004, was consistent with that 
derived from the model including parallax and or- 
bital motion, -0.0153±0.0004. 

4.9. Second Order Effects 

With the slope parameter included we had ac- 
counted for all the physical effects which we ex- 
pected to be present in the lightcurve. Having 
found that the residuals showed no further varia- 
tion at a level detectable above the photometric 
noise, we did not attempt to include second order 
effects such as xallarap etc. 

4.10. Final Model 

All our models are compared in Table [2] and 
the parameters of the best-fitting models are pre- 
sented in Table [3] and Figure Q] We plot the x 2 
for each link in the MCMC chain for all parame- 
ters against one another in the best-fitting model 
in Figure [5j This plot was used as a diagnostic 
throughtout the fitting process, as any correlations 
between parameters display distinct trends as the 
chain moves towards the minimum. The caustic 
structure changed during the course of event, so 
in Figure |6] we show the structure at two distinct 
times; the first at the time of the first caustic cross- 
ing during the anomaly, and the second at the time 
of closest approach. 
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Table 2 

The x 2 for the best-fitting model in each class, comparing 
m > and u < solutions. 



Model 


x 2 






wo > 


u < 


Standard 


6331.950 


6331.950 


Standard+Parallax 


6064.773 


6099.487 


Standard+Parallax+Orbital Motion 


5544.071 


5690.003 


Standard+Parallax+Orbital Motion with Slope 


4782.367 


4802.606 


Standard+Slope 


5334.378 


5334.378 



Table 3 
The best fitting parameters. 



Parameter 
(Units) 


Standard 


Standard 
+Slopc 


Standard 
-(-Parallax 


Standard+Parallax 
+Orbital Motion 


Standard 4-Parallax 
+Orbital Motion+Slope 


x 2 


6331.950 


5334.378 


6064.773 




5544.071 


4782.367 


Ax 2 " 


1549.583 


552.011 


1282.406 




761.704 


0.0 


to (HJD') 


5344.32 


5344.38 


5344.47 




5344.83 


5344.69 


±0.01 


±0.01 


±0.02 




±0.03 


±0.02 


M 


0.4089 


0.403 


0.403 




0.381 


0.386 




±0.0009 


±0.001 


±0.001 




±0.001 


±0.001 


*E (d) 


43.82 


44.84 


43.49 




43.4 


44.3 


±0.08 


±0.09 


±0.09 




±0.1 


±0.1 


so 


0.7692 


0.7725 


0.7717 




0.7792 


0.7750 




±0.0005 


±0.0005 


±0.0006 




±0.0007 


±0.0007 


1 


0.0705 


0.0677 


0.0695 




0.0683 


0.0654 




±0.0005 


±0.0005 


±0.0006 




±0.0006 


±0.0006 


«0 


0.180 


0.171 


0.198 




0.297 


0.221 




±0.003 


±0.003 


±0.003 




±0.006 


± 0.007 


p 


0.01963 


0.01912 


0.01931 




0.0163 


0.0165 




±0.00006 


±0.00007 


±0.00008 




±0.0001 


±0.0001 








0.18 




0.96 


0.37 








±0.02 




±0.04 


±0.05 








-0.124 




0.09 


0.01 






±0.007 




±0.01 


±0.01 


ds/dt (yr _1 ) 










0.53 


0.49 












±0.02 


±0.02 


da/dt (yr _1 ) 










-1.21 


-0.37 












±0.06 


±0.08 


Slope (mag/yr) 




-0.0160 








-0.0153 






±0.0004 








±0.0004 


a All timcstamps 
b Improvement in 


are abbreviated to H,ID'= HJD-2450000.0. 
X 2 relative to that of the best fitting model. 
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Fig. 5. — Chi squared contours plotted as a function of the parameters fitted in the MCMC fit for the best 
model. The red, orange, green, light blue, dark blue and purple colors indicate the regions with Ax 2 <l-6cr 
(respectively) from the best-fit solution. Inset: Close-up of the coutours for the parallax parameters. 
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5. Physical Parameters 

The purpose of this model is to ultimately ar- 
rive at the physical parameters of the lens and 
source, which can be achieved using the known re- 
lations between these and the lensing parameters 
obtained from the modeling. 

Chiefly of interest is the mass of the lensing sys- 
tem, Mf^tot, which can be determined explicitly 
for events where parallax is measurable provided 
the angular extent of the Einstein radius, #e, is 
known from: 



c 2 AU 6 E 
AG tte' 



(7) 



where is the Einstein radius projected from the 
source onto the observer's plane. The model pa- 
rameter p represents the angular size of the source 
star 9s in units of the angular Einstein radius 9e- 
We derive this from the crossing time taken for the 
source to travel behind the lens, = fi6s where 
p is the relative source-lens proper motion, p can 
then be written as: 



(8) 



= e_s 

These parameters also yield the distance to the 
lens, 



AU 



7TL = #ette 



which in turn yields the projected separation be- 
tween the lens components: 



aj_ = soD L t 



(10) 



and the relative proper motion between lens and 
source, when combined with t^: 



/' 



(11) 



The appreciable lens orbital motion during this 
event also allows us to test whether the companion 
object is bound to the primary lensing mass, via 
the ratio of its kinetic to potential energy: 



KE/PE 



(s Re) 3 7 2 
8tt 2 Ml 



(12) 



where 7 relates the two lens orbital parameters: 
-f 2 —(ds/dt/so) 2 + (da/dt) 2 and where the masses 
are in units of Mq, the distances in AU and time 
measured in years. 

However, these expressions include two key terms 
which are as yet, unknown: 9s, the angular source 
radius and Ds, the distance to the source. In order 
to extract the physical characteristics of the lens, 
we therefore turned our attention to the charac- 
teristics of the source. 

5.1. Source Star 

Long-exposure V, I images were acquired by 
the CTIO 1.3 m at several epochs which enabled 
us to plot the color-magnitude diagram for the 
field including the source star (Fig. [7]). By observ- 
ing the event at different levels of lensing magni- 
fication, these data can be incorporated into the 
model which yields the source and blended light 
fluxes, f s , fb for those data, and hence the instru- 
mental magnitudes and colors of the source and 
blend. But we note that these uncalibrated fluxes 
also suffer from the high degree of extinction along 
the line of sight to the Galactic Bulge. To calculate 
the dereddened color, (V — I)s,o, and magnitude, 
Is,o, of the source, we needed to calibrate the in- 
strumental fluxes f s and fb relative to a standard 
candle. 

Fig. [7] clearly shows a locus of stars centered 
at /Rc,inst=15.821±0.05mag, (V - J)RC,mst=- 
0.350±0.05mag. This consists of a clump of red gi- 
ant stars, for which stellar theory predicts a stable 
absolute luminosity, varying only slightly with age 
and chemical composition. Their frequent occur- 
rence^ ^inakes_these objects useful as standard can- 



dles. IStanek et al.l ([19981 ) established photometric 



calibrations for red clump magnit udes which were 
later refined bv lAlves et ail (l2002h using Hipparcos 
data. Most recently. iNataf et al.l ( 2012 ) were able 
to measure the dereddened apparent magnitude 
of the red clump stars at the Galactocentric dis- 
tance, /rc,o= 14.443. By mapping the distances, 
Drc, to red clump stars in the Galactic Bar as 
a function of Galactic longitude, I they found an 
apparent viewing angle on the Bar of 0Bar=4O°, 



i?o sin <p + I 



Drc 



= cos I + sin I cot (f>, (13) 



where Nataf et al. (|2012l) measured i?o to be 
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8.20 kpc. For the field of MOA-2010-BLG-073 
(Galactic coordinates: Z=4.81030°, 6=-3.50131°), 
we derive .0^(7=7.48 kpc, and we assumed that 
the source star lies behind the same amount of 
dust as the Red Clump stars and at the same 
distance. Scaling the dereddened apparent mag- 
nitude of the red clump stars, /rc,0j appropriately 
for the slightly closer distance of the stars in this 
field, /Rc,ap P = Irc,q+ A/, where: 



A/ = 51og 10 R /D 



RC- 



(14) 



We found Al = 0.20 mag, and so the distance 
modulus to the red clump a nd the source i n this 
field is /i?c,app= 14.24 mag. iBensbv et~ai1 (|201ll ) 
determined the intrinsic (V — I)Rc,o=l-09mag for 
red clump stars, so we were able to derive their ab- 
solute V magnitude of Mv,rc,O=0-97 mag. Com- 
bining these results with the measured A(V— T)i ns t 
and AVinst between the source and red clump in 
the CTIO data, we then derived the dereddened 
color, (V — I)s,o and magnitude Is,o of the source, 
summarized in Tabled 

Bessell k Brett! (|l988h provided a relationship 
between (V — I) and (V — K) color indices, and 
Kervella et all (|2004[) related (V - K) to angu- 



lar radius for giant and dwarf stars. Having thus 
determined Dg and 9s, we derived the physi- 
cal parameters of the lensing system, which are 
also summarized in Table |U The color and large 
source radius of 14.7±1.3i? Q implies this star 
is a K-type giant, which is consi stent with the 
observ ed photometric variability. I Jorissen et al 
(J1997J) found that red giant stars with spec- 



tral types later than early-K are all variable, 
with amplitudes increasing from microvariabil- 
ity to several magnitudes towards cooler tem- 
peratures and t imescales from days to years. 
Kiss et all ( 20061 ) notes that irregular photomet- 
ric variability may be caused by large convection 
cells, or may actually be the result of a number 
of simultaneous periodic pulsation modes, and 
many examples have been identi f ied from time- 
doma i n surveys ( Wrav et al. 20041 Wozniak et al 



2004t lEver k Blakei 120051 : lciechanowska et al 
20061 e.g.). We note that the star was de- 



tected by 2 MASS (jSkrutskie et al.l 120061 ) as source 
2MASS J18101138-2631226 with colors (J - 
#)=0.76±0.078mag and (iJ-if s )=0.284±0.079mag 
Although the 2MASS field is crowded the star's 



PSF is distinct and their photometry for it has 
the best-quality AAA flag. These colors are 
consistent with a giant star and with our de- 
rived value for 1^0=12.345 mag, when we take 
i nto account AATg=0. 24mag from the vvv survey 
(iGonzalez et alj|2012h . 



Nataf et al.l (|2012l ) explain that their value of 



our viewing angle of the Galactic Bar is a "soft up- 
per bound" because the distance along the plane 
to the greatest density of stars along the line-of- 
sight to a triaxial Bar structure is less than the 
distance to the structure's major axis on the far 
side and greater on the near side. As the physical 
parameters derived for MOA-2012-BLG-073 arc 
somewhat dependent on the v alue of our yiewin; 
angle of the Galactic Bar and iNataf et al.l (|201 
emote consistent results with values as low as 25°, 
we explored the potential impact of this on our 
results. A reduced viewing angle would produce a 
smaller distance to the source, changing it's dered- 
dened magnitude and color. The resulting increase 
in source radius produces a corresponding reduc- 
tion in the value of Dz,and increases in the lens 
masses. However we found that the physical pa- 
rameter values do not change by more than the 
errors quoted in Table 21 implying that this is not 
the dominant source of uncertainty. Finally, we 
computed the physical parameters derived from 
the best-fitting uq < Omodel for comparison and 
found that the masses derived changed by < la. 

6. Discussion 

The Earth's movement during this relatively 
long time scale (tE=44.3d) microlensing event re- 
sulted in a gradual shift in our perspective on 
the lensing system, breaking the symmetry of the 
caustic. Meanwhile, the change in projected sep- 
aration of the lensing objects modified the shape 
of the caustic just as the source's trajectory hap- 
pened to pass close by. If not for these subtle vari- 
ations, it can be seen from Figure [5] that a source 
trajectory uq > would produce exactly the same 
light curve as a uq < trajectory. As it is, for 
this event, the uq > solution best explains our 
observations, though the difference in \ 2 relative 
to the corresponding uq < solution is very small 
(Ax 2 =20.2) compared with the \ 2 of both fits. 

The measurable parallax signature enables us 
to determine the masses of the lensing bodies. The 
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Table 4 

The physical parameters of the lensing 
system and source star, derived from 
the model including slope, parallax and 
orbital motion, plus color information. 



Parameter 


Units 


Value 


e s 


/nas 


9.143±0.792 




mas 


0.557±0.09 


Rs 


Rq 


14.7±1.3 


M L ,i 


Mq 


0.16±0.03 


M La 


M } 


11.0±2.0 


M L , tot 


Mq 


0.17±0.03 


D L 


kpc 


2.8±0.4 


a± 


AU 


1.21±0.16 


KE/PE 




0.079 


Proper motion 


mas yr _1 


4.60±0.4 


-Zs.inst 


mag 


15.554±0.007 




mag 


15.335±0.007 


(V - 7) S , inst 


mag 


-0.22±0.01 


^RC.inst 


mag 


15.821±0.05 


(V — J)RC,inst 


mag 


-0.350±0.05 


Irc,o 


mag 


14.443 


(V - 7)r C ,o 


mag 


1.09 


Is,o 


mag 


13.976 


V s ,o 


mag 


15.197 


(V - I) s ,o 


mag 


1.221±0.051 


(V-K) s ,o 


mag 


2.852 


K s ,o 


mag 


12.345 


J (2mass) 


mag 


13.686±0.053 


H (2MASS) 


mag 


12.926±0.057 


Ks (2mass) 


mag 


12.642±0.054 
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primary lensing object has Mi i i=0.16±0.03M Q , 
making it an M-dwarf star. The companion's 
mass is A/ L ,2=H.0±2.0Mj. This places it the 
brown dwarf desert, though we note that this 
traditionally refers to close-in companions, and 
since microlensing and direct imaging measure 
only the projected separation, we know only 
their minimum orbital semi-major radii. Regard- 
less, it is clear that MOA-2010-BLG-073Lb is 
close to the mass threshold for deuterium burn- 
ing (~0.012M Q =12.6Mj) quoted as the nominal 
boundary between pla nets and brown dwar fs es- 
tablished by the IAU (jChabrier et al.l l2005h . So 
what kind of object is it? 

No further orbital or metallicity information is 
available for this event, which might have shed 
light on its evolutionary history. Theoretical 
isochrones predict that a star of this low mass will 
not have lost a significant amount of material over 
its lifetime, so we can say that MOA-2010-BLG- 
073L b formed as a high mass-ratio binary. Models 
of protoplanetary disks have the expectation that 
disk mass, -Mp, will scale line arly with star mass 
M* (IWilliams fc Cieza 11201 lh . M D /M* - 1% at 
young ages, but this would limit Mjj - and hence 
M P , max - to -1-10 Mj in the case of MOA-2010- 
BLG-073L. So it seems questionable whether such 
a massive companion could have formed in a pro- 
tostellar disk, via either core accretion or gravita- 
tional instability. 



Bonnell et al.l (|2008D ; lKroupa fc Bouvierl (|2003h 



discuss a number of mechanisms which can pro- 
duce an M-dwarf/Brown Dwarf binary following 
gravitational fragmentation in a molecular cloud: 

1. Embryo rejection model the nascent binary 
was ejected from a dynamically unstable 
multiple protostellar system, leading to the 
loss of its accretion envelope before the sec- 
ondary component could acquire enough 
mass to become a star. 

2. Collision model the binary was prematurely 
ejected from a larger protostellar system by 
the close passage of another star. 

3. Photo-evaporation model the accretion enve- 
lope around the binary was photo-evaporated 
by the nearby presence of a massive star in 
the birth cluster before the secondary could 



accrete enough mass to become a star. 

4. "Star-like" model the object formed as a 
normal stellar binary with low-mass compo- 
nents. 



The embryo rejection model predicts (jBate. Bonnell fc Bromm 



2002) that the maximum separation of binaries 
surviving this process is a max ~ 10 — 20 AU. We 
cannot rule out this scenario as we only measure 
the projected separation of the lens, which is nev- 
ertheless below a max . Since MOA-2010-BLG-073L 
is a field object we have no information regard- 
ing the proximity of other stars during its birth, 
so the collisional and photo-evapouration mod- 
els are equally plausible. However, we note that 



Whitworth fc Stamatellosl (|2006l ) derived a mini- 
mum mass for primary fragmentation components 
of 0.004M Q ~4.2Mj. As this threshold is below 
the mass of MOA-2010-BLG-073Lb the simplest 
explanation is that the companion is an extremely 
low-mass product of star formation. However, 
we note that the mass ratio of this system would 
match that of an object in the brown dwarf desert 
if the host were a mo r e mas si ve star. Recent sur- 
veys (jDieterich et al.l (|2012f ); lEvans et all (|2012 , 
e.g.), though restricted to massive companions 
with a > 10 AU) hint that the brown dwarf desert 
may extend to M-dwarfs and beyond 3 AU, which 
would make MOA-2010-BLG-073L a rare example 
of its type. 
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Fig. 4. — The OGLE-II photometry was binned 
using a series of bins of varying width in time. 
The data in each bin was fitted with two functions 
a) constant brightness b) with a linear slope, and 
the RMS around these curves was averaged over 
all the bins in the lightcurve. The average and 
maximum RMS for each binning is plotted here 
against the width of the bin, indicating variability 
over different timescales. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Fig. 6. — Maps of the caustic structure during 
the anomaly and at the time of closest approach. 
The red line and arrow indicate the trajectory of 
the source in a reference frame centered on the 
barycenter of the lensing system, while the red 
dots indicate the position of the source at these 
times. The green and blue dots indicate the po- 
sitions of Mi ! (largest) and M^ 2 at both times 
(radii of dots not to scale). 
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Fig. 7. — The instrumental color-magnitude di- 
agram for the field of view including the lensing 
source star. The position of the source is marked 
with a blue star, relative to the center of the Red 
Giant Clump highlighted with a red diamond. The 
green circle indicates the color and V magnitude 
of light blended in the CTIO photometry. 
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Fig. 8. — The ratio of planet mass to host star 
mass plotted against semi-major axis for all exo- 
planets for which these parameters are a vailab le 
(source: exoplanet.eu Schneider et al. ( 201 lh ) . 



The true orbital semi-major axis is plotted where 
it has been measured, otherwise the projected 
(minimum) separation has been indicated with an 
arrow, where the base of the arrow marks the mea- 
sured value. Very few objects have been found be- 
tween log 10 (<?)=- 1 - -2. This corresponds to the 
brown dwarf desert. The location of MOA-2010- 
BLG-073L, highlighted in cyan and bold, lies be- 
tween the brown dwarf and planetary regimes. 



